Rientrant metallicity in the Hubbard model: the case of honeycomb nanoribbons 
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Based on the Cluster Perturbation solution of the Hubbard hamiltonian for a 2-D honeycomb 
lattice we present quasiparticle band structures of nanoribbons at half filling as a function of the 
on-site electron-electron repulsion. We show that at moderate values of e-e interaction ribbons 
with armachair shaped edges exhibit an unexpected semimetallic behavior, recovering the original 
insulating character only at larger values of U. 
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The repulsive interaction among electrons is respon- 
sible of the failure of single particle picture and of the 
opening/widening of energy gaps in solids. The Hubbard 
model is the paradigm to describe this phenomenon: suf- 
ficiently large values of the on-site e-e repulsion inhibit 
the inter-site hopping favoring in this way an insulat- 
ing behaviour. The 1-atom thick 2D honeycomb lattice 
(graphene) does not contradict this picture: many body 
effects due to on-site Coulomb repulsion have been shown 
to lead, for sufficiently strong interactions, to semimetal- 
to-insulator transition 1 * 2 - as well as to other deviations 
from Fermi-liquid behavior such as unconventional quasi- 
particle lifetimes^, long range antiferromagnetic order— 
and spin liquid phased. 

In this paper we show that for honeycomb nanorib- 
bons the repulsive e-e interaction may be responsible of 
a metallic phase in ribbons that in the single particle 
picture are semiconducting. This appears to be another 
extraordinary property of the 2-D honeycomb lattice. 

It is well known that honeycomb nanoribbons mani- 
fest peculiar properties related to the topology of their 
edges^: according to single-particle theory ribbons with 
armchair shaped edges may exhibit a finite energy gap 
depending on their width 7 - - — , while ribbons with zigzag 
edges are metallic and become insulating only after the 
inclusion of an antiferromagnetic orde r 10 ! 11 . The modi- 
fications of the single particle band structure of zigzag 
graphene ribbons due to e-c interaction has been in- 
vestigated within a mean field solution of the Hubbard 
mode l 12 i 13 , showing spin polarization of edge states and 
gap opening at the Fermi level. How the single parti- 
cle band picture evolves to a quasi-particle one and how 
the electronic states of both armchair and zigzag honey- 
comb ribbons are modified by on-site Coulomb repulsion 
described as a true many body term is the question we 
address in this paper. 

We have adopted a many body approach based on the 
Cluster Perturbation Theoryii (CPT). CPT belongs to 
the class of Quantum Cluster theories^ that solve the 
problem of many interacting electrons in an extended 
lattice by a divide- and- conquer strategy, namely solv- 
ing first the many body problem in a subsystem of finite 
size and then embedding it within the infinite medium. 



Quantum Cluster theories represent some of the most 
powerful tools for the numerical investigation of strongly 
correlated many-body systems. They include Dynami- 
cal Cluster Approach 1 ^, Cellular Dynamical Mean Field 
Theoryi 7 . as well as CPT and have found an unified lan- 
guage within the variational scheme^ based on the the 
Self Energy Functional approach-^. CPT has many inter- 
esting characteristics and gives access to non trivial many 
body effects in a relatively simple way: it exactly repro- 
duces the limits U/t = (non- interacting band limit), 
U/t = oo (atomic limit) and the Mott-Hubbard metal- 
to-insulator transition for intermediate values of U/t; it 
recovers most of the characteristics of the exact solution 
of the 1-dimensional cas o 14 ' 20 ; finally it is relatively easy 
to implement and, at least for simplest systems, without 
much numerical effort. In its variational for m 18 ' 21 (Vari- 
ational Cluster Approximation (VCA) ) it can be applied 
to systems with spontaneously broken symmetry describ- 
ing antiferromagnetism 2 - 2 - and superconductivit y 23 ! 24 . In 
zigzag honeycomb ribbons VCA has been applied to 
study the transition from topological to antiferromag- 
netic insulator—. 

In CPT— the lattice is seen as the periodic repetition of 
identical clusters (Fig. [1] ) and the Hubbard Hamiltonian 
is written as the sum of two terms, an intra-cluster (H c ) 
and an inter-cluster one (V) 

6 = Hi + Hu> = H c + V 

I Ijtl' 

where the summations are over all the clusters and 

ijcr i 
ijcr 

Since in the Hubbard model the e-e Coulomb interac- 
tion is on-site, the inter-cluster hamiltonian V is single 
particle and the many body term is present in the intra- 
cluster hamiltonian H c only, a key feature for the prac- 
tical implementation of the method. Having partitioned 
the Hamiltonian in this way we may write the resolvent 

operator G as G^ 1 = z — H c — V = G c — V and from 



this 

G = G c + G C VG (1) 
The one-particle propagator 

S(knw) = < ^olC^kr^o > 

+ < *o|cknG4nl*0 > 

is obtained exploiting the transformation from localized 
to Bloch basis 

4„ = ^E<( k )*^ k - (R!+n) 4 

and similarly for Ck n - Here a™(k) are the eigenstate coef- 
ficients obtained by a band calculation for a superlattice 
of L identical clusters, identified by the lattice vectors Ri, 
each cluster containing M sites at positions n; n is the 
band index and the summation is over the N = L X M 
lattice sites. We get 

Giknc) = ^^ e - k -^-^')|ar(k)| 2 fe(M (2) 

ii' 

where Gw (kw) is the superlattice Green function, namely 
the Fourier transform of the Green function in local basis 

IV 

This is the quantity that can be calculated by cq[T] 

Q iV (koj) = g% v [lu) + J2 Bij (kw)S,v (kw) (3) 

3 

where M x M matrix Hy(ka;) is the Fourier transform of 
G C V involving neighboring sites that belong to different 
clusters. Once the cluster Green function in the local 
basis Qly (lu) has been obtained by exact diagonalization, 
eq. [3] is solved quite simply by a M x M matrix inversion 
at each k and to. The quasi particle spectrum is then 
obtained in terms of spectral function A(k.Lu) 

ACku) = - Y~ ImGfknuj). 

n 

The key approximation in this derivation has been to 
identify the many electron ground states of the extended 
lattice as the product of cluster few electron ones. This 
is certainly wrong except at U/t = or U/t = oo and at 
intermediate values of U jt it is important to verify the 
accuracy of the results by using larger and larger cluster 
sizes. In practice this procedure is limited by the dimen- 
sions of Hilbert space used in the exact diagonalization, 
dimensions that grow exponentially with the number of 
sites. 

Fig. Q] illustrates the ribbon geometries and the open 
chains that are periodically repeated to reproduce both 
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FIG. 1: (Color online) Geometric configurations of honey- 
comb ribbons with zigzag and armchair edges identified ac- 
cording to the notation of reference^. The single (a,c,d,e) and 
double (b) chains used to tile the extended ID lattice are also 
shown. 

zigzag and armchair ribbons. We have considered zigzag 
ribbons of different widths and in the case of ribbon 
< 1.1 > both single chain (4 sites) and double chain 
(8 sites) tilings in order to test the influence of cluster 
size on the quasiparticle spectrum. In Fig. [5] we show 
the QP Density Of States (DOS) obtained for a specific 
value of U jt and the two above mentioned cluster sizes, 
compared with the corresponding non-interacting DOS. 
We notice that the main features (peak positions, gap 
opening) do not depend much on the cluster size and we 
may confidently use the smallest cluster size in all cases. 
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FIG. 2: (Color online) Density of Quasi Particle states ob- 
tained for U/t = 4 for a zigzag terminated ribbon assum- 
ing different cluster sizes (Fig [1] a,b) compared with non- 
interacting DOS. 

As a result of the inclusion of e-e interaction extra 
structures appear in the quasi particle spectrum be- 
low (above) filled (empty); this is what happens also 
in real materials where Hubbard correlation may be re- 
sponsible of severe energy renormalization, quasiparticle 
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FIG. 3: (Color online) Value of the energy gap as a function 
of U /t showing a drastically different behavior in zigzag and 
armchair ribbons. For zigzag termination results obtained for 
different cluster sizes (single and double chain, Fig. Q](a) and 
(b) ) and for different ribbon widths are shown. 
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quenchin g 26 i 27 and the appearance of short lived satel- 
lites structures 2 ^. The main effect of Hubbard correla- 
tion however is the opening, for sufficiently strong inter- 
action U > U c , of a well defined gap. This is due to the 
well known Mott-Hubbard mechanism: U inhibits double 
occupancies of sites and in this way makes electron hop- 
ping from site to site less and less energetically favorable, 
driving the system across a metal-to-insulator transition. 
This is what happens in model systems at half occupation 
and in real materials 2 ^. As shown in Fig. [3] the value of 
U c for zigzag ribbons depends on the ribbon width. Pre- 
vious VCA calculations^ have found for wider ribbons 
U c « 3t, a results not too far from the present one taking 
into account the larger ribbon width and the variational 
procedure used in that work. Recent ab-initio estimate of 
the screened on-site Coulomb interaction in graphene 3 ^ 
provide values of the the same order (U/t ~ 3.5). 

For the armchair terminated ribbon the situation is 
completely different. As previously mentioned, according 
to single particle theory ribbons with armchair shaped 
edges may be metallic or insulating depending on their 
width£~— . In particular the ribbon of interest here ex- 
hibits a finite gap in the non-interacting picture and we 
may expect the one-site e-e interaction to reinforce this 
insulating behavior enlarging the gap. This, surprisingly 
enough, is not quite the case: switching on e-e correlation 
the gap first diminishes, reaches zero at U/t = 3 and only 
after that grows linearly with U/t (Fig. [3]). At U = U c 
the QP band dispersion (Fig. 2]) becomes linear around 
k = and the system semimetallic. 

We may look at the character of quasiparticle states 
by considering the local spectral function, namely 

A(k,i,ui) = i£„ /m £( knw ) I <( k ) I 2 - We notice that 
the effect of e-e repulsion is to delocalize states close to 
the gap region that for U/t = are localized either at the 
edges or in the inner part of the ribbon. In particular for 
U/t — 3 the localization of gapless states responsible of 
the semimetallic behavior is equally distributed across 
the ribbon. We have found that it is possible to rcpro- 



FIG. 4: (Color online) Quasiparticle band structure for arm- 
chair terminated ribbon <1.1> obtained with U/t = (a),= 2 
(b),= 3 (c), = 4(d). The black line superimposed in (c) is the 
band structure obtained in the single particle picture assum- 
ing different hopping terms (see text). 



duce the same gapless band dispersion with the same de- 
localized character in a single particle picture attributing 
different hopping terms to edge sites and to inner ones: 
this is shown in Fig. EJc) assuming hopping between edge 
states exactly twice as big as hopping between inner ones. 
In this sense the net effect of the on-site e-e repulsion is to 
renormalize the inter-site hopping selectively across the 
ribbon and to make hopping between doubly coordinated 
sites (edge sites) more favorable than hopping between 
sites with triple coordination (inner sites) . The gap clos- 
ing would then be due to a mechanism similar to a me- 
chanical strain. The analogy between mechanical strain, 
e-e repulsion and magnetic field as responsible of gap 
tuning in graphene has been recently investigate d 3 ^ 32 
showing in particular that in the mean field approxima- 
tion a local nearest neighbor Coulomb interaction may 
create in a 2D honeycomb lattice non trivial magnetic 
configurations and metallic phases with broken time re- 
versal symmetry 3 ^. In the present study we have shown 
that also on-site Coulomb repulsion if treated as a full 
many body term may induce a semimetallic behavior in 
a semiconducting f D honeycomb lattice. 

What is surprising is that the Hubbard interaction has 
opposite effects in the two types of honeycomb ribbons: 
it opens a gap in the metallic ribbon and closes it in 
the semiconducting one. A similar behavior has been 
reported in a recent paper on correlation effects in topo- 
logical insulators 2 ^. Starting from an extended 2D hon- 
eycomb lattice made semiconducting by spin-orbit (SO) 
interaction the authors show that for U/t = 3 the existing 
gap closes down. We have repeated the calculation using 
the same SO parameter and calculated also in this case 
the gap as a function of U/t. The results are shown in 
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FIG. 5: (Color online) Energy gap as a function of U/t for 
2D honeycomb lattice with (circles) and without spin-orbit 
(diamonds) interaction. 



Fig. [5] for a 2D honeycomb lattice with and without SO 
and again in the semiconducting system we find a regime 
of U/t where the energy separation between filled and 
empty states decreases and the pre-existing gap closes 



down to zero. And this even if the original semiconduct- 
ing behavior has a completely different physical origin 
(SO interaction in 2D instead of armchair termination in 
ID). 

In conclusion, we have found that many body effects 
associated to local e-e repulsion may be responsible of a 
(semi)metallic phase in systems with honeycomb lattice. 
The Hubbard mechanism that inhibits double occupan- 
cies of sites, instead of reducing the ability of electrons to 
jump from site to site, induces a selective renormalization 
of inter-site hopping and the energy separation between 
filled and empty states becomes zero at a specific k-point. 
This behavior appears to be a characteristic of honey- 
comb topology - or perhaps more generally of bipartite 
lattices - at half occupation and suggests the existence of 
an [/-dependent additional symmetry^. The Electrons 
trapped in artificial lattices with honeycomb geometry^ 
are the best candidates where this anomalous behavior 
induced by e-e repulsion (gap closing/opening for semi- 
conducting/metallic systems) can be experimentally ver- 
ified. 
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